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An exact solution representing black holes in an expanding universe is found. The black holes 
are maximally charged and the universe is expanding with arbitrary equation of state {P = wp 
with — 1 < yw < 1). It is an exact solution of the Einstein-scalar-Maxwell system, in which we 
have two Maxwell- type U(l) fields coupled to the scalar field. The potential of the scalar field is an 
exponential. We find a regular horizon, which depends on one parameter (the ratio of the energy 
density of U(l) fields to that of the scalar field). The horizon is static because of the balance on the 
horizon between gravitational attractive force and U(l) repulsive force acting on the scalar field. 
We also calculate the black hole temperature. 

PACS numbers: 



Black holes play a central role in astrophysics 
as well as string theory, and may possibly be cre- 
ated at the LHC [H . Not surprisingly, they have 
been studied intensively over the past 35 years. Im- 
portant progress has been made in understanding 
the Hawking process [2] and the states responsible 
for black hole entropy at the microscopic level Q . 
However many problems remain unresolved: does 
cosmic censorship [4] hold?, what happens when 
black holes collide?, how does accretion of mat- 
ter affect the thermodynamics of black holes? and 
how does it affect the growth of black holes? To 
answer these dynamical questions one needs time- 
dependent solutions of the Einstein equations con- 
taining black holes. In this letter we shall focus on 
solutions representing black holes in a background 
FLRW universe. 

There have been many previous attempts to 
obtain black holes embedded in the Friedmann- 
Lemaitre- Robertson- Walker (FLRW) universe. 
The Einstein- Straus model is perhaps the simplest 
one It is a patchwork of Schwarzschild black 
holes with an FLRW universe. However these 
black- holes are time symmetric, and so they do not 
describe a dynamical black hole in the Universe. 

One well known black hole candidate in FLRW 
universe is the McVittie solution P, which is a 
spherically symmetric, time-dependent solution of 
the Einstein equations. The solution approaches an 
FLRW universe at "infinity," and looks like a black 
hole near the "horizon". However, as shown in 
the McVittie solution is disqualified as a black hole 
(or a point mass singularity) in the FLRW uni- 
verse. Recently, Sultana and Dyer constructed a 
more sophisticated black hole solution in a dynam- 
ical background by a conformal technique [8] . The 
matter content is null dust and ordinary dust. The 



solution tends to an Einstein-de Sitter spacetime 
asymptotically. This model, however, violates the 
dominant energy conditions. 

Assuming self-similarity, we can show that a reg- 
ular black hole may exist only in an accelerating 
universe, but this requires numerical study p. The 
analytic solution found by Carr and Hawking de- 
scribes a self-similar spacetime with a regular black 
hole but it approaches asymptotically a "quasi" 
FLRW spacetime which has a deficit angle, but not 
an exact flat FLRW spacetime [l^. There are also 
discussions of "dark energy" accretion onto a black 
hole in a universe. 
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If a cosmological constant is present, we have 
the Schwarzschild-de Sitter (SdS) and Reissner- 
Nordstrom-de Sitter (RNdS) solutions ^,d3i]. Al- 
though these spacetimes are static, they may be 
converted by a coordinate transformation into the 
form of a black hole in an exponentially expand- 
ing Universe [1^. Multi-black hole solutions in 
a de Sitter Universe were found by use of ex- 
tremely charged black holes and their collision dis- 
cussed [m, [l6|. This Kastor-Traschen (KT) so- 
lution is a time-dependent generalization of the 
Majumdar-Papapetrou solution, which describes 
extremely charged RN black holes [l^ • Similar so- 
lutions were given in [Til, Ell . 

Another time-dependent cosmological black hole 
system was found from the compactification of in- 
tersecting brane solution in higher dimensional uni- 
fied theory [i^. As clarified in [2l| the global pic- 
ture of dynamical solution describes a multi-black 
hole system in the expanding Universe filled by 
"stiff matter" . We shall call it the MOU solution. 

Here we generalize these two solutions and 
present an exact solution describing a cosmological 
multi-black hole system with an arbitrary power 
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law expansion. This is a solution of general rela- 
tivity with a scalar field with an exponential po- 
tential and two Maxwell- type U(l) fields coupled 
to the scalar field. The solution has regular event 
horizons, approaches asymptotically an exact flat 
FLRW spacetime without a deficit angle, and no 
singularity exists outside the horizons. 

The above known solutions take the following 
form: 
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The KT solution with TV-black holes located at the 
coordinate position r(^) = 1, • • • is given by 
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where (5(i) is the charge of the i-th black hole, and 
the scale factor of the background Universe is given 
by a(f) (X exp(i^of) {Hq: constant). It is asymp- 
totically de Sitter spacetime. The MOU solution 
discussed in [20,, 21l given by 
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where Qr^i) and Qs{i) are the conserved charges of 
time-dependent and static branes, respectively, and 
the scale factor is given by a = (f/fo)^^^, which also 
holds for an expanding universe with stiff matter. 

By changing the time coordinate, these solutions 
can be rewritten in the form of a brane system, 
discussed in detail in [20], as 



where 
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Here ut and ns are appropriate no n- negative real 
numbers with the constraint + = 4, and to is 
a constant. The transformation of the time coordi- 
nate is given by t/to = a^/^^(f), where to is fixed 
as to = ^o"^ foi" the KT solution and 3fo/4 for the 
MOU solution, respectively. Setting = 4, we 
find the KT solution, while the case with = 1 
corresponds to the MOU solution. 



Assuming the metric form (|3]) ~ (|6]), if we take 
an arbitrary real value for (or n^), we find that 
the scale factor a in the form of ([1]) is given by any 
power function, i.e., a (xt^^ where p = n^^/n^. (We 
regard p = oo as an exponential expansion.) 

When the universe expands with an arbitrary 
power by a scalar field, one needs an exponential- 
type potential [22|. In fact, if we have the universe 
filled by a scalar field with the potential 



V Voex.p{-atv^) , 
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the scale factor increases as a oc f ^ with p = 2/a^, 
where Hi^ = SttGn is the gravitational constant. 
Thus we should choose a = = \/2nslnT- 

We shall therefore adopt the following action: 

^ = ^d^x^g\^^R-\g^^(d^^){d,^) 
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where g^^y is a spacetime metric, $ is a scalar 
field with the potential V{^) given by (O, and 
fI^J{I = S,T) are two Maxwell-type U(l) fields, 
which couple to the scalar field with the coupling 
constants A/. The vector potentials are given by 
A\i\ and n/ are their degeneracy factors. 
The metric (|3j) with (jH) plus 



^^ = -^mns/2 \n{HT/Hs) , 
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with (O, (O and a^^Fq^o = ^t(^t - 1)/4 is really 
an exact solution of the system ([8|), if we assume 



At = \/2ns/nT , \s = -\/2n^r 



(11) 



and + = 4 [23]. For = 4 and = 1, we 
recover the KT and MOU solutions, respectively. 

The above solution with arbitrary gives a 
multi-black hole system in an expanding universe 
for which the scale factor and effective equation 
of state are given by a ex with p = n^^/ng, 
and P = wp with w = |^ — 1, respectively. 
Note that w takes an arbitrary value in the range 
of —1 < w < 1, corresponding to the value of 
1 < < 4. 

We summarize some typical solutions in Table |T1 
In order to discuss the spacetime found here in 
detail, in what follow, we consider a single black 
hole system. For simplicity, we assume that two 
charges are equal, i.e., Qt = Qs =• Q- We shall 
rewrite the metric (|3]) as 



ds' 



-r'^U-^dP + {dr^ + r'^dftl) (12) 
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with dH) and Ht = t ^ Hs = l^r 

= to/Q, ds'^ = ds^ IQ^, r = r/Q, and i = t/to are 
dimensionless variables. The metric ([12]) depends 
on only one parameter r, whose physical mean- 
ing is given as follows: The energy density of the 
scalar field is uniform at t = to, which is given 
by p^{to) = 3n^/16to. While the total density of 
the U(l) fields is evaluated on the horizon for the 



static extreme RN black hole with the charge Q 
as Pu(i)|// = For the time-dependent black 

hole, p^\h and Pu(i)|/f are different from the above 
values, but their orders of magnitude are still the 
same. Thus we can claim that r is related to 
the ratio of two energy densities at the horizon as 



~ Pv(i)/p^\h (see 21, 2^). The limit of r ^ oo 
gives the static extreme RN black hole. 
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TABLE I: Some values of the typical parameters of a black hole system in a universe, for which expansion law and 
the equation of state are given by a oc t ^ and P = wp, respectively. 



The circumference radius R = QR^ which is a 
geometrically invariant variable, is given by 

R = rU = (l + tf)^ + . (13) 

The curvature singularity appears at f = — 1 and 
f = — where R vanishes. Analyzing the be- 
haviour of trapping horizons in the limit of f ^ 
and near horizon geometry as in [21], we find the 
horizon radius or R-) which satisfies the fol- 
lowing equation: 

r{R^ -l) = Rl, T{Ry -1) = -Rl. (14) 

The spacetimes are classified by their causal 
structure into three types: Type I < 2), Type 
II {ut = 2), and Type III (n^ > 2). 

In Type I, there are two horizons, and 
which are the roots of Eqs. (p!4|) . Since r > 0, 
we find R-^ > 1 > R- > 0. We show the horizon 
radii in terms of r in Fig. [T] (a). For Type II, if 
r > 1 there are two horizons, = \Jt jir — 1) 
and R- = \/t/{t ^ 1), but if r < 1, we find only 
one horizon, R-. In Type III, if r > Tcr, we find 
two roots ^+,1 and ^+,2 (^+,1 < ^+,2) for the 
equation for Rj^ (p!4j). where Tcr = n^^^'^{nT — 
2)-(nT-2)/2^2. From the detail analysis of space- 
time structure, we find there exist two horizons; R- 
and = ^+,1 {R- < 1 < ^h), but = ^+,2 
is not a cosmological horizon except for the KT 
solution [2^. The cosmological horizon turns out 



to be time-dependent just as that in an accelerat- 
ing universe [25|. On the other hand, if r < Tcr, 
we find one horizon, R-. If r = Tcr, there are two 
horizons, R- and = i^+,2 (degenerate) (see 

Fig. [11(b)). As we see from Table IH Type I and 
III correspond to a decelerating and accelerating 
universes, respectively. 




(a) riT = 1 (b) TT-T = 3 



FIG. 1: Two horizon radii {R- and or Rn) for 
riT = 1 and 3. i?H and i?c, which is not a cosmological 
horizon, degenerate at Tcr = 3v^/2. 

Since the present spacetime is spherically sym- 
metric and the near horizon is "static", we can 
calculate the surface gravity (see [21] for details). 
Hence we find the black hole temperatures on the 
horizons (Tbh) by the surface gravity n± as 

We depict the behaviour of the temperatures in 
Fig. [21 They are finite and vanish in the limit of 
r cxD, i.e., the extreme RN spacetime. 
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and the Japan-U.K. Joint Research Project. 



(a) Ut 



(b) Ut 



FIG. 2: Black hole temperatures on two horizons (R- 
and or Ru) for tit = 1 and 3. Tbh = QTbh- 

It may be interesting to discuss the thermody- 
namics because we can define the entropy and tem- 
perature in these time-dependent spacetimes. We 
can easily extend the present solution to arbitrary 
dimensions |24|. The details including the analysis 
of global structure and study of thermodynamical 
properties will be given elsewhere [24]. Some ques- 
tions for future work are : 

(1) Can we find more realistic black hole solutions? 
It may be straightforward to include rotation (cf 
[Toj ). This is under investigation. As for neutral 
black holes in the Universe, it may be difficult to 
obtain the analytic solution because such a system 
is non-BPS state even in the static case, and the 
radius of a black hole increases in time due to ac- 
cretion of matter . 

(2) Can we extend black hole thermodynamics to 
time-dependent spacetimes? The present time de- 
pendent solution may provide a good tool for ana- 
lyzing this question. 

(3) Can we discuss some dynamical process with 
the present or extended solutions? Black hole col- 
lision can be discuss in the contracting Universe 
(to < 0) just as the KT spacetime. We can also dis- 
cuss the brane collisions with multi time-dependent 
branes, which is generalization of (26| . 

(4) Some solutions have a link to intersecting brane 
systems in higher-dimensional supergravity model. 
If Ut is a non- negative integer, we may regard 
Ut and Us as numbers of branes. It is true for 
Ut = 1, in which case we can derive the four- 
dimensional effective action ([8]) from compactifica- 
tion of the time-dependent M2-M2-M5-M5 brane 
system in 11-dimensional supergravity theory (see 
Appendix in |2l|). Hence it may be interesting to 
see whether there is any fundamental or deep rea- 
son for this link. 

Work along these lines is in progress. 
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